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1 Introduction 

Research interest in controllable two level systems, which have been enthusi- 
astically called quantum bits or qubits, has grown enormously during the last 
decade. Behind a huge burst of activity in this field stands an idea of what 
is possible in principle but extremely difficult to achieve instrumentally - the 
fascinating idea of quantum computing. The very principle of quantum super- 
position allows many operations to be performed on a quantum computer in 
parallel, while an ordinary 'classical' computer, however fast, can only handle 
one operation at a time. The enthusiasm is not held back by the fact that ex- 
ploiting quantum parallelism is by no means straightforward, and there exist 
only a few algorithms (e.g., [1, 2]) for which the quantum computer (if ever 
built) would offer an essential improvement in comparison with its 'classical' 
counterpart. Even if other uses of quantum computing prove limited (which 
might or might not be the case), its existence would most certainly lead to a 
breakthrough in simulations of real physical many-particle systems. 

Whether or not the ultimate goal of building a working quantum computer 
is ever achieved, both experimental and theoretical studies of properties of sin- 
gle or entangled qubits are flourishing. One of the most compelling reasons for 
this is an exciting overlap of the possibility of a future technological break- 
through and the reality of research in fundamentals of quantum mechanics. 
There exist various experimental realisations of the qubit, amongst which 
solid state qubits are of particular interest as they provide one of the most 
promising routes to implementing a scalable set of qubits, which is one of the 
minimal requirements for implementing quantum computations. However, any 
solid state qubit, albeit representing effectively a two-level system, comprises 
a huge number of internal degrees of freedom whose unavoidable coupling to 
the environment leads to loss of coherence. Quantum computations require a 
set of fully or at least partially entangled states on which some unitary opera- 
tions are performed. Decoherence would make evolution of states non-unitary 
and would lead to the loss of entanglement between the states. Thus, the 
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loss of coherence before a sufficient amount of quantum operations was per- 
formed would be the major impediment in using solid-state qubits in quantum 
computations, ft is believed that tens of thousands of unitary operations are 
required for quantum computation to become a reality [3] so that sufficiently 
long decoherence times (much longer than those currently archived in the best 
solid state qubits) should be achieved experimentally. It necessitates a better 
theoretical understanding of realistic mechanisms of decoherence. 

In these lectures, after illustrating in a simple way the main features of 
decoherence in a generic qubit coupled to the environment (Section 2), we will 
consider onset of decoherence in one particular realisation of qubit, namely 
charge Josephson junction (JJ) qubit. First, we will describe briefly what is 
the JJ qubit (Section 3) and then focus on a mechanism widely believed to 
be responsible for one of the main channels of decoherence for the charge 
JJ qubit, namely its inevitable coupling to fluctuating background charges 
(Section 4). This model has been thoroughly investigated in all regimes [4-7] 
and has a tutorial advantage of being exactly solvable within a fully quantum 
approach [7] and showing a rich variety of different regimes with a non-trivial 
dependence on temperature and on (unfortunately ill known experimentally) 
the strength of coupling between the qubit and the fluctuating charges. 

In Section 5 we will offer a solution for decoherence rate in this model 
which is formally exact at an arbitrary temperature T. In the 'high-T regime' 
(which could still correspond to rather low temperatures), the decoherence 
rate saturates and becomes T-independent, while at low temperature it turns 
out to be linear in T and behave non-monotonically as a function of the 
coupling strength between qubit and the environment. In conclusion we will 
also consider the relaxation rate, albeit only perturbatively with respect to the 
coupling strength, and demonstrate that the model can qualitatively explain 
the experimentally observed [8] quasi-linear behaviour of the spectral density 
of noise with humps at certain frequencies. 

2 Coupling to the environment and decoherence 

Before considering a realistic solid-state qubit, we start with illustrating what 
is the loss of coherence in a generic qubit. Such a qubit is a two-level system 
so that its Hamiltonian can be mapped to that of spin \ and written 



where B is an effective 'magnetic field' (measured here in energy units). States 
of the qubit can be described in terms of its density matrix, 
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When the qubit (or any system) is in a pure state, one can always find a basis 
where p = \i) In any other (rotated) basis, the density matrix of a pure 
state obeys p 2 = p. In a rotated basis, the density matrix performs a unitary 
evolution described by the Heisenberg equation of motion, 

(3) 

whose solution is p{t) = U p{0)W where in this trivial case the evolution oper- 
ator U = c~ lHot . In a semiclassical language, such an evolution is convenient 



at 



H ,p(t) 




to visualise using the Bloch sphere representation, shown in Fig.l. There the 
spin evolution is parameterised by two angles on a unit sphere, 6 and 4>. For 
the simple system described by Eqs. (1) - (3) above, such an evolution is 
simply rotating around z axis, provided that the 'field' B = (B z ,0), and the 
initial state of the spin was diagonal in a different basis. The Bloch angles are 
related to the density matrix by 

p(t) = i(l + S(t).<x). (4) 

An inevitable coupling to the environment (bath) can be schematically 
described by the Hamiltonian 

H = H o + Koupi + iibath ■ (5) 

Now it is the full density matrix of the qubit + bath which obeys the Heisen- 
berg equation of motion. As we are interested in states measured on the qubit, 
we need to introduce a so called 'reduced' density matrix, p q = Trbath P where 
the trace is taken over all the bath states with the proper Gibbs weight. We 
will describe a consistent way of performing such a trace in Section 5. Here we 
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want just to illustrate how the coupling to the environment makes the time 
evolution non-unitary and leads to the loss of coherence. 

Let us consider a model in which there is a minimal, so called longitudinal 
coupling of the qubit to the bath, 

Koupi = Xa z . (6) 

Let us further assume that the qubit was prepared in a pure state but in a 
different basis, so that its density matrix has both diagonal and off-diagonal 
elements, 

fpn PiaWn / 
VP21 P22 J \f 1 - n 

with |/| 2 = n — n 2 att = assuming that the qubit was prepared in a pure 
state, p 2 = p. Let finally switch off the x component of the 'magnetic field' 
B. In the absence of coupling, the spin that represents our qubit on the Bloch 
sphere of Fig. 1 would simply rotate at a constant frequency u>o oc B z around 
z axis. The longitudinal coupling effectively introduces and additional, fluc- 
tuating time dependent component of the field, B z {t). This fluctuating field 
results from the thermal noise in the bath and inevitably destroys coherence 
as we will now show. 

The spin rotation is described in a scmiclassical language by the Landau- 
Lifshitz equation, 

dS „ „ 
— = S x B. 

At 

In the absence of the transverse component of the magnetic field, the diagonal 
elements of the density matrix, related to S by Eq. (4), remain constant. In 
the presence of the environment-induced time-dependent component of the 
longitudinal magnetic field, S rotates with a changing frequency resulting in 
the off-diagonal elements of p acquiring the following time-dependence: 

Pl2(t) = \ (S x +iS y ) = p 12 (0)e^ t+i /o ^(r)dr (7) 

The contributions due to the fluctuating field, B z (t), should be averaged 
over time. It is natural to assume that fluctuations of B z are Gaussian, as 
they are due to an very large number of degrees of freedom in the bath. For 
Gaussian fluctuations, one can use the standard averaging formula, 

(e^)=e-^ 2 >, 

so that 

(pi2(*)>=Pi2(0)e iwo * (e^o^Wj 

■l*dTj*dT'(B(r)B(r')) 



= Pi2(0)e iwot exp 
= Pl2 (0)e^ t exp 



, J* dr(B{Q)B( T ] 



= pi 2 (0)e^ ot e- r2t (8) 



Decoherence and relaxation in charge qubits 



5 



Here we have semi-formally defined the decoherence rate B 2 as the rate of 
relaxation of the off-diagonal part of the density matrix, using that the corre- 
lation function of the fluctuating field B z depends only on the time difference. 
In this definition, F 2 can still be time-dependent. However, for large enough 
t, longer than the longest relaxation time for thermal noise in the bath which 
makes B z fluctuating, F 2 should saturate at some limiting value. Thus, making 
a bit more formal definition, we arrive at 

r 2 = ±- = lim / (It(b(O)B(t)) = S B (w « 0) , (9) 

where 5s (w) is the noise power spectrum of the fluctuating field B, 

S B (co)= [ (B{0)B(t))e- iujt dt. (10) 

J — oo * ' 

The above calculation should be considered just as an illustration as we have 
just introduced Gaussian fluctuations of the 'magnetic field' B z . However, it 
is very easy both to make this sort of calculations fully quantum-mechanical, 
and to generalise the model beyond the longitudinal coupling. 

To show that the effective coupling to the bath leads to the thermal noise 
that directly results in the appearance of decoherence given by Eq. (9), one 
needs to model the bath. The most standard theoretical approach to such 
modeling and thus to decoherence by the environment is based on spin-boson 
models [16, 17], where the environment is modelled as a set of harmonic os- 
cillators, 

£bath = 5> k &]A (11) 
k 

The linear coupling in Eq. (6 ) should be understood as coupling to all the 
oscillator degrees of freedom, 

Koupi = Xa z = (\bl + h.c.) a z (12) 

k 

The fluctuation-dissipation theorem allows one to express the noise power for 
the coupling operator X(t) (in the interaction picture) as 

Sx(uj) = (Xl) + (Xi u ) = 2JMcothg, (13) 

where J{lo) depends on the density of oscillator states, i.e. on the spectrum 
Ct>k in Eq. (11), and on the coupling A (ui) in Eq. (12). The decoherence rate 
is still given by Eq. (9), where Sx{ui) should be substituted for Sb(w). For 
the most typical, Ohmic model, J (ui) awso that _T 2 oc T. As we will show 
later, such a linear dependence on temperature is characteristic for a realistic 
model to be considered later in these lectures, but only for sufficiently low T. 
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It is straightforward to generalise our considerations beyond the longitu- 
dinal model. To this end one should take into account the cr^-proportional 
contribution in the qubit action - this would be general enough even if one 
leaves the coupling as it stands in Eqs. (6) or (12). 

3 Charge Josephson junction qubit 

Not any two-level system could serve as a qubit. There is a set of requirements 
that such a system must satisfy in the first place. In 1997 David DiVincenzo 
[9] formulated five criteria which have to be satisfied by a physical system 
considered as a candidate for quantum computation, conditions which are 
widely known as the 'DiVincenzo checklist': 

i) One needs well-defined two-state quantum systems (qubits); 

ii) One should be able to prepare the initial state of the qubits with sufficient 
accuracy; 

iii) A long phase time coherence is needed, sufficient to allow for a large 
number (~ I0 4 ) of coherent manipulations; 

iv) Control over the qubit 's Hamiltonian is required to perform the necessary 
unitary transformations; 

v) A quantum measurement is needed to read out the quantum information. 
At the present time any known physical system falls short of the lsit stipulated 
above. Currently lots of different possibilities of building a quantum computer 
arc investigated including nuclear spins, quantum dots, Josephson junctions, 
trapped ions, optical lattices, electrons on liquid helium and some others with 
each of them being quite far from satisfying this or that condition of the above 
set. 

A route which seems to be one of the most promising is to build a system 
of Josephson junction (JJ) qubits: the qubits which are based on Josephson 
junctions and utilise the charge and flux degrees of freedom. The main ad- 
vantages of the superconducting devices involving Josephson junctions are a) 
a relative easiness to manufacture - the lithographic methods used to fabri- 
cate them are well established; b) controllability - gate voltages controlling 
charge Josephson qubits can be adjusted with a very high degree of accuracy 
allowing, in particular, the high-fidelity preparation of the initial state of the 
qubit; and finally c) measurability - the techniques for almost noninvasive 
measurements which can be used for this type of qubits are quite advanced. 
It is point iii) in the DiVincenzo checklist - decohcrence - which causes the 
greatest worry for the J J qubit, as well as for any other realisation of a solid 
state qubit. We stress again that 'spin up' and 'spin down' states in solid state 
qubits are formed by a substantial number of electrons (which is not the case 
for nuclear spins which can be effectively decoupled from the external world) . 
All degrees of freedom interact with the environment thus causing the loss of 
quantum coherence. Whether this problem of decoherence can or cannot be 
satisfactorily solved is most likely to determine the future of JJ qubits. 



Decoherence and relaxation in charge qubits 



7 



The JJ qubit is made of superconducting islands separated by a Joseph- 
son junction; the scheme is depicted in Fig. 2. In the scheme, 'SC marks the 
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Fig. 2. Scheme for charge JJ qubit 



superconducting islands and the hatched area is the Josephson junction with 
capacity Cj and Josephson coupling Ej. The superconductors are separated 
from the leads by capacitors C, which do not allow any tunneling, and these 
capacitors are biased by some controllable gate voltage V g . Cooper pairs can 
tunnel through the Josephson junction thus changing the total charge of each 
island. 

The charge JJ qubit is a two-level system which states are different by 
the charge of a single Cooper pair. The equilibrium number of Cooper pair in 
each island is controlled by the gate voltage. When it is tuned to be close to 
half integer ~ uq + 1/2 , the states with uq and no + 1 Cooper pairs (circled 
in Fig. 3) can be arbitrarily close to each other in energy while the distance 
to other energy levels (crossed in Fig. 3) is of the order of the charging energy 
that - in the temperature units - could be tens or even a hundred degrees 
of Kelvin which is much higher than typical experimental temperatures. We 
assume that the Josephson coupling energy is much smaller than the charging 
energy so that the states are really discriminated by charge. (In the opposite 
limit, a so called flux J J qubit can be - very successfully - built). Thus we 
effectively separate all the charge states with n < no or n > no + 1 and thus 
have a two level system, whose states are made of a macroscopic number of 
electrons. In future we will be referring to these as states 'spin up' and 'spin 
down'. 

More formally, the total energy of the system which consists of charging 
and Josephson contributions can be written as follows: 

s -2c>Tc( 2| ^-^) 2 -^ cos0 ' (14) 

where & is the superconducting order-parameter phase shift between the is- 
lands. In properly chosen units, the pair of n and & are canonically conjugated 
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Charging energy 




Fig. 3. Two charge states 



coordinate and momentum since they generate correct Hamiltonian equations 
of motion. Differentiating with respect to momentum dn/dt — dE/dO = 
Ej sin produces the correct equation for the Josephson current (in the units 
chosen Ej coincides with the critical current), and differentiation with respect 
to coordinate d0/dt = -dE/dn = (CV g - 2n)/{2Cj + C) = Vj produces the 
correct value for the voltage drop Vj across the junction. The quantisation of 
(14) then gives: 



H = 



1 



2C., + C 



CV a 



Ej COS O, 



,0 



(15) 



After taking the projection of the Hamiltonian (15) onto two states of 
interest, in this 'spin up' and 'spin down' basis the charging part can be 
written as (SE c /2)a z , where SE C is the charging energy difference between 
the two states: 



SE C (Vg) 



2C., + C 



n 



CV n 



(16) 



The commutator between n and dictates the commutation relations [n, e ±ie ] 
-p e ±i6> f rom w hlch it immediately follows that in the 'spin up' and 'spin down' 
basis cos is given by a x /2. Summarizing the above we substitute the Hamil- 
tonian (15) with its projection: 



. 5E c {V g ) 
tl = cr. 



Ej 



(17) 



As the result we have a controllable two level system or a qubit. The control 
is exercised by changing the gate voltage V g on which the a z coefficient SE C 
depends (16). If 6E C is maintained large 5E C ^> Ej, which is characteristic 
for charge J J qubit, then the Josephson part of the Hamiltonian is irrelevant 
and evolution amounts to acquiring phase shift between 'spin up' and 'spin 
down' states in the initial mixture. On the contrary, tuning charging part 
to the degeneracy point SE C = stimulates the spin flip process. It can be 
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shown that the ability to switch for arbitrary time to the above two regimes 
is enough in order to provide all necessary one qubit quantum operations or 
quantum gates. 

In building a quantum computer, one and two-qubit gates will be required. 
The details of how to build them are beyond the scope of these lectures and 
can be found in the comprehensive review [10] of Makhlin et. al. But what 
is crucial for any quantum computation is to maintain entanglement between 
the qubits. This will be undermined by inevitable decoherence. Therefore, 
understanding its mechanisms and studying how it sets in is crucial for future 
progress in this area. 

4 Decoherence - Fluctuating Background Charges model 
4.1 The model 

The main problem with charge JJ qubits (as with any other solid state qubits) 
is that they lose quantum coherence too quickly due to unwanted, but un- 
avoidable, coupling to the environment. Recent experiments [ll-15]show that 
hundreds of elementary quantum operations can be achieved before coher- 
ence is destroyed. This is at least two orders of magnitude short [3] of a rough 
estimate of what is required for non-trivial quantum calculations. Thus the 
problem of decoherence is currently the major obstacle for the progress of 
qubits of this type, and requires a close theoretical examination. 

Spin-boson models will not offer us any insight into how real physical pa- 
rameters of the problem affect decoherence. To obtain this we need a physical 
model describing the real processes which decohere the qubit. It is widely be- 
lieved [11, 12] that in charge J J qubits the main contribution to decoherence 
comes from coupling of the qubit to some charge fluctuators present in the 
environment. The fluctuating background charges (FBC) model, suggested in 
this context in [4] and then also used in [5, 6, 7, 18], is an attempt to build a 
microscopic model for such a dominant channel of decoherence. 

The schematic picture of how the qubit interacts with the FBC is shown in 
Fig. 4. Impurities which contribute to decoherence sit on the substrate which 
the nanocircuit was grown on, close enough to a metallic lead to allow tun- 
neling to and from the lead. Solid circles in Fig. 4 are charged impurities and 
transparent circles are their image charges. These dipoles interact electrically 
with the qubit, which in turn behaves like a dipole, since one of the super- 
conducting islands has some number of excessive Cooper pairs and the other 
one lacks the same number. The interaction depends on the number of such 
extra pairs and thus is different for the two charge states of the qubit, enforc- 
ing dependence of the coupling term on a z . Although a static dipole-dipolc 
interaction could only shift the qubit states without any loss of coherence, 
the tunneling between the impurities and the metallic lead makes charges on 
the impurities to fluctuate, effectively creating a random time-dependent field 



10 A. Grishin, I. V. Yurkevich, I. V. Lerner 



substrate 



localised state 



metallic gate 




Fig. 4. Interaction between the qubit and fluctuating background charges 

on the qubit. This random field causes decoherence. This is similar to the 
generic case considered in Section 2: the time-dependent electric field coupled 
to a z plays the role of the random 'magnetic field' B z in Eq. (8). In the fol- 
lowing section, we will show that within the Fermi golden rule approximation 
(valid for a weak coupling between the qubit and the FBC), both decoher- 
ence and relaxation can be calculated in a way similar to that schematically 
developed in Section 2, with an additional advantage of calculating the noise 
power-spectrum within a fully microscopical model. In this section, we restrict 
considerations to the longitudinal model of pure decoherence. Its advantage 
is that one can go well beyond the Fermi golden rule, and calculate the deco- 
herence rate (but not the relaxation!) non-perturbatively and within a fully 
quantum description valid for any temperature, including T — > (in which 
limit there is no decoherence in the absence of relaxation). 

The Hamiltonian corresponding to the described model can be written as 

H= 6 -^a z -^-a x +a z V + H B ; V = ± £ v t d\d t ; 

i 

Hb=J2 £ HH + E M U + h.c] + E • (18) 

i i.k k 

Here d f , d\ are the operators of annihilation and creation of an electron on 
the i-th impurity; c k , c k are the operators of the conduction electrons in the 
metal; iki are the hybridisation amplitudes; e° is the energy of the localised 
state on the i-th impurity; £k are the energies of the conducting electrons; Vi 
is the coupling strength between the qubit and the i th impurity. 

The charge J J qubit corresponds to 8E C » Ej, which is always the case 
for the small enough capacities (i.e. small superconducting islands and the 
JJ junction in Fig. 2. A non-trivial circuit is supposed to consist of many 
qubits, and most of the time each particular qubit is in an idle regime for 
which the condition above is fulfilled. During relatively short times necessary 
for operations involving spin- flips the longitudinal model, Ej = 0, is not 
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appropriate. However, if there is any hope for a working qubit, the onset of 
decoherence should only happen in the idle regime. 

Apart from temperature, there are three other parameters of dimensions 
of energy for each fluctuator: the coupling strength v i} tunneling rate 7, (7* = 
2tt ^ k |£ki| 2 <5(w— £k)), and fluctuator energy (note that the bare energy e° is 
renormalised by hybridisation). All three parameters are broadly distributed, 
and temperature can be considered 'low' for some of them and 'high' for 
others. In the Section 5 we will show that these two regimes are defined as 
follows: 



T < min 
T ^> max 



or T <C 7i low-T regime 



and T ^ 7^ high-T regime 



(19) 



In the high-temperature regime, the decoherence rate can be calculated 
classically [4, 5] which gives for one fluctuators (omitting the index i) 



r 2 = 4 



7/2, w>7 
w 2 /47 , v < 7 



(20) 



The high temperature treatment of the model would be justified if the high 
energy impurities, \si\ > T, were frozen. However, the same hybridisation 7 
which makes impurity charges to fluctuate, broadens their energy positions. 
Such a broadening creates a Lorentzian tail which is, of course, power-law 
suppressed but nevertheless allows for a contribution from the energetically 
remote impurities. Moreover, even impurities with \si\ <C T may not sat- 
isfy the high-T inequality (19) and contribute in a non-classical way. We will 
present below some arguments that the impurities most relevant for decoher- 
ence are likely to to be in the low-temperature regime. 



4.2 Why low temperature? 

Let us first illustrate why the hybridisation leads to the Lorentzian tail. To 
simplify notation, we will be considering a single impurity. The quadratic 
Hamiltonian of the bath Hb can be diagonalised by means of some lin- 
ear transformation [19], after which it becomes Hb = En e n^i"n • The in- 
teraction part V gets transformed into a sum over all exact states: V = 



(«/2)£n, 



u u ct ] a 



, where coefficients u n are related to the density of 



states (DoS) v e {u>) on the impurity broadened by hybridisation with the con- 



duction band: 



v e {£n)5; here 5 is the level spacing in the conduction 



band, and the DoS is given by the Lorentzian 



7/2 



7t (u - e) 2 +7 2 /4 ' 



(21) 
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Fig. 5. Qubit is coupled to the transitions between all exact states of the bath 
Hamiltonian Hb, with weight being proportional to u* n u m . Coefficients Un,u m are 
governed by the Lorentzian |u„| 2 = v e (e n )8. 

In this notation, low temperature manifests itself in smallness of the inter- 
action term V , to which only n and m in the vicinity of the Fermi level will 
effectively contribute (see Fig. 5). Since the smallness of u coefficients is con- 
trolled by the Lorentzian (21) (rather than distribution function), the effect of 
such a high energy impurity can only be power law suppressed. 

The second observation concerns very stringent geometrical restrictions on 
where an 'effective' impurity can be and a broad distribution of energies of 
impurities. In order to allow tunneling to and from the impurity, it should not 
be too far away from the metallic gate, which implies that the total volume 
available for 'effective' impurities cannot exceed Ifim x (lOAi?) 2 , where 1/xm 
is a characteristic size of a metallic gate and Af is the Fermi wavelength in 
the metal. If we are hoping to find at least one low energy \ei\ <T impurity in 
this volume, the total number of impurities there should be of order of D/T, 
where D is the width of the distribution function of energies. For temperature 
which varies in experiments between 30mK and 50mK, and for D having a 
typical chemical value of leV the above ratio is of order of 10 5 . Having that 
many impurities in a volume which in total accommodates ~ 10 6 atoms is 
completely unrealistic. 

So, the contribution of high energy impurities is suppressed but not 
strongly enough in order to exclude them from consideration, and low energy 
impurities effectively do not exist. If for a particular sample there happened 
to be one, this sample will exhibit orders of magnitude stronger decoherence, 
and should better be discarded. 

The condition |e,| » T enforces low temperature regime which requires a 
full quantum mechanical treatment. Such a treatment, valid at any tempera- 
ture, will be offered in the next Section. 
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5 Exact solution for FBC model at arbitrary 
temperature 

5.1 Calculations 

To calculate the decoherence rate we will need to know how the off-diagonal 
elements of the reduced density matrix of the qubit decay with time. The 
dynamics of the full density matrix of the whole system qubit + bath is given 
by the Hciscnbcrg equation of motion: 



dpit) _ ■ 
dt 



H,p(t) , (22) 



where H is the Hamiltonian given by the Eq. (18) with Josephson energy 
being taken equal to zero. Since in the longitudinal (Ej = 0) model the 
interaction term commutes with the qubit Hamiltonian, one can trace out 
the bath degrees of freedom thus yielding the following result for the reduced 
density matrix of the qubit: 



P#(0) p$(0)e-^D(t) 
p^(0)e^D*(t) p$(0) 



P (q) (t)= m^^'^'m^ ^' I , (23) 



where the standard separable initial condition for the full density matrix, 

p(0) = p iq) {0) ® pb\ Pb = e~ /3J * B / Tre ~ /3 * B > was assumed. The diagonal 
elements of the reduced density matrix do not evolve since there are no spin flip 
processes present, and the dynamics of the off-diagonal elements is governed 
by the decoherence function D{t): 

D{t) = ( e i(H B +v)t e -i(H B -v)t\ ) ( 24 ) 



where averaging should be performed with the bath Hamiltonian H B . At long 
times t the decoherence function must decay exponentially, D(t) ~ c~ r2 *, so 
that the decoherence rate r 2 is defined as 

T 2 = -3?c lim t- 1 In D(t) , (25) 

in the agreement with Eq. (9). 

The average in Eq. (24) can be represented as the following functional 
integral with the Grassmann fields £ and r/ defined on the Keldysh contour: 



D(t) = Z- 1 / VCT>£Vr)*Vrie e e ^ * 



s [(] iSoM 2/ I>< (*')£<* (*')£< (t')dt' 



X 



-i / E[t kl <(*')C I (*')+*k l C(*')%(*')]dt' 

x e " K k ' ; . (26) 
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Fig. 6. Dependence of Vi(t) on time along the Keldysh contour 



Here Z is the same functional integral but with Uj(t') = 0, the fields £ s , £* 
correspond to the localised state on the i-th impurity, ?7 k , ?y£ to the conduction 
electrons, the impurity 5o[£] an d rree electrons So[rj] actions are given by: 

So W = / J2vUt%d t , - e k ) Vk (t')dt' , (27) 

OK k 

and the Keldysh time dependent coupling Vi(t') is zero everywhere on the 
contour apart from the interval (0, t) where on the upper branch it takes the 
value of Vi, and on the lower branch it is equal to — Vi (see Fig. 6). Since the 
action So[r)} is quadratic and the hybridisation term is linear in r\ fields, it is 
straightforward to integrate out the conduction electrons: 



D(t) = Zq 1 j veve, 



iS [«] S «r(*l)^i(*l,*2)? 3 (t 2 )dtldt 2 

e e " K 13 



i / I>i(ti)£i*(*i)U*i)dti 
e <=* 1 (28) 



Again, here Z is the same functional integral but with Vi(t') = 0. All the 
information about the conduction electrons is now encoded in the self-energy 
matrix which is defined on the contour: 



(29) 



with gk(ti,t2) being the Green function of free electrons. The full action in 
(28) is quadratic and therefore the corresponding functional integral can be 
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written in a standard symbolic trace-log notation: 



D(t) 



,Trln 1+GC/ 



(30) 



where both G and U are matrices in the space of impurity indices, and both 
depend on two times along the Kcldysh contour. The symbol Tr here stands 
for trace with respect to the impurity indices and for integration along the 
Kcldysh contour over all the times involved. The Green function G obeys the 
following integro-differential equation: 



G m j(t' ,t 2 ) — 8ij6(tl, t 2 ) ■ 



(31) 



with 5-function 8{t\, t 2 ) being defined on the contour. The impurity matrix U 
is made out of the time-dependent coupling Vi(t'): 



Uij(ti,t 2 ) = 5ij8{t u t 2 ) 



Vi{t 2 



(32) 



Being written explicitly, the symbolic expression in (30) becomes 



D(t) = e 



- E 

n = l 



(-IT 



tr / ■■• / dtidt' 1 ...dt n dt' n G(t 1 ,t' 1 )U(t' 1 ,t2)---G(t n ,t' n )U(t' n ,t 1 ) 

2„ integrals (33) 



where tr stands for a trace over impurity indices only. 

Since the matrix U{ti,t 2 ) as function of its both times is non- vanishing 
only between and t, all the contour integrals in (33) can be converted into 
ordinary integrals from to t by introducing the Keldysh structure for both 
matrices G and U. Both arguments of G(ti,t 2 ) can either be on the upper 
or on the lower branch of the Keldysh contour leaving four choices which 
generate the following Keldysh structure: 



Gij(ti,t 2 ) tl 2ecK i ^ Gij(h — t 2 )t ly2 e(-oc,+oo) 1 



K g>(ti-t 2 ) g ij {t 1 -t 2 ) / 



(34) 



where Qij (t\ — t 2 ) and Qij (t\ — t 2 ) arc the time-ordered and anti-time-ordered 
Green functions. Using the Keldysh structure for the [/-matrix, 

Uij{ti,t 2 ) tl 2ecK i ^ \8ij a z S(ti - t 2 ) Vi(t 2 ) ti ^g^^^ooj , 

and the general rule, 



I A(t u t')B(t',t 2 )dt' i ^ J i(ii,t')^5(t',i 2 )di' 



ti,26( — oo,+oo) 
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we rewrite the expression (33) for the decoherence function as follows: 

- E ^S-S dt 1 dt 2 ..At n g(t 1 -t 2 )(v®i)g(t 2 -t 3 )(v®i)...g(t n -t 1 )(v®i) 
D(t) = e " =1 '° 

(35) 

Now the symbol Tr stands for the trace over both impurity and Keldysh 
matrix indices, while (v <g> 1) denotes the matrix with elements Vi on the main 
diagonal in the impurity space and which is the unity matrix in the Keldysh 
space. 

It is convenient to do a standard rotation [20] in Keldysh space: 

Qij = La z QijD; L = -^=(1 - ia y ) , (36) 

which converts the four-element matrix in (34) into the three-element one of 
retarded, advanced and Keldysh Green functions: 



In the rotated basis the expression (35) for the decoherence function becomes: 
D(t) = e 



- E t-Sr- Tr/---/ dtidt 2 ...dt„(«lgl* x )S(ti-t2)(fi®*x)S(t2-i3)...(«®*x)S(tn-tl) 



(38) 

Since due to the time translation invariance each G depends only on the differ- 
ence of its time arguments, the n th order integrand depends onn-1 differences 
in times, while the integration over the last time variable produces the overall 
factor of t. The region of integration becomes, at arbitrary time, quite com- 
plicated. However, when time t is much bigger than the characteristic time on 
which Q{t) function decays (which is true if t ^> ,T _1 ), all the integrals 
over the time differences can be extended to the entire axis. Then the integral 
has a convolution structure in time and, performing a Fourier transform and 
a straightforward summation that restores the logarithm, it finally reduces to 

* / ^Trl„[i+I(«®^)a(^)] 

D(t) = e -~ . (39) 

Using the definition of decoherence rate (25) and calculating trace in the 
Keldysh space explicitly we have: 



/J 
— trln 
2tt 



i + -G^)--G»-G» 



(40) 
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where tr refers only to the impurity matrix indices. 

It should be stressed that the analytical structure of the expression in the 
r.h.s. of Eq. (40) ensures that there is no decoherence at zero temperature. 
Indeed, using the standard relation [20] between the components of Keldysh 
Green functions at equilibrium, G K (ui) = tanh(/3w/2)[£r (w) — G A {w)\, which 
at T = reduces to G k (uj) = sgn(w)[G fl (w) - G a (uj)}, we have: 



r 2 (T = o) 



+ OC 

-Jfte / -^trln 

ZTT 



(41) 



Separating the retarded and advanced parts in the expression above, and then 
expanding the logarithm in both of them, we see that starting from the second 
order all the integrals are convergent at u) — > oo and therefore equal to zero 
since the retarded (advanced) Green function is analytic in the upper (lower) 
half plane. Combining the first order contributions from both parts we see 
that it is purely imaginary and therefore also vanishes. 

Assuming that the impurity states are well localised, the tunneling ampli- 
tudes iki can be written as tki = fjVol _1 / 2 e lkri , where is the position of 
the i-th impurity. This in turn means that the self energy matrix (29) has the 
following representation: 



0- 



(42) 



Since the characteristic distance / between the impurities is much bigger than 
the Fermi wavelength, I » Xp , all the off-diagonal elements of the self energy 
matrix will be averaged out while diagonal can be represented as 



£ JwHiiTrVltwI^a; ) 

^— ' Id — £k 



L k 



c*i(w) 7i(w) 



±i- 



(43) 



This allows one to solve (the Keldysh rotated analogue of) Eq. (31): 



- zf R H 



U) - £. 



lu - ei =F i7»/ 2 



(44) 



where both the real and imaginary parts of the self energy matrix were as- 
sumed kj-indcpendent, and the real part was absorbed by renormalisation of 
the energy levels £j = e^ + ai/2. In such a diagonal approximation for the self- 
energy, the contributions from all the impurities are independent and the total 
decoherence rate is given by their sum. Substituting Eq. (44) for the Green 
functions into Eq. (40), subtracting (to improve the integral convergency) the 
identically zero expression (41) for F and taking the real part of the resulting 
expression, we obtain the following contribution of a single impurity at energy 
Sj = e to the decoherence rate: 
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+ OC 

.<*■> - - / : 



du> 
4ir 



In 



1 - 



w 2 7 2 cos\i- z (lu/2T) 



w 2 7 2 + 4[( w _ e )2 + 7 2/ 4 _ w 2 / 4 ]< 



(45) 



5.2 Discussion of the results 

It is convenient to introduce an auxiliary 'spectral function' A(iv), 

yl ( w ) = T—u ^7 \2i 2 2 ' e± = £± ^V 7 ^ 2 - 7 2 , (46) 

1 + 4(w — £+) (w — £-Y/v z ^ z 2 

thus rewriting Eq. (45) as follows: 



duj 
Air 



In 



1 — A(iv) cosh 



2T 



(47) 



Behaviour of /l(o>) is qualitatively different depending on whether g = 
v/'y < 1 or g > 1. In the former case vl(w) is peaked at w = e, while in 
the latter the peak splits into two Lorentzians which are centered at uo = e± 
and travel in different directions away from the ui — e point as v increases. 
This signals a qualitative change between the weak coupling (<? <C 1) and the 
strong coupling (g » 1) regimes. 

The decoherence rate is determined by the overlap of the temperature 
function, cosh _2 (o;/2T), and A(u>) in (47). The former is centered at u = 
and exponentially decays for > T, while the latter has the power-law decay 
away from the above described peaks. If the temperature function is wide and 
the Lorentzian peaks are sitting inside T, the overlap is solely controlled by 
A(u>) as if the temperature function was equal to 1 everywhere. In the oppo- 
site case when the temperature function is narrow either with respect to the 
distance to the yl-function peaks or to their width, the overlap is controlled 
by cosh~ 2 (w/2T) and A(u>) can be replaced by A(u> = 0). We will be call- 
ing these two cases high-temperature (T ^> max|e±|,7) and low-temperature 
(T <C min|e±| or T <C 7) regimes correspondingly. 

In the high-tcmpcraturc regime, T 3> max|£±| , 7, the decoherence rate 
saturates [4, 5] at e-independent value: 



r (hi g h) = _ I ^ ln[1 _ m] = 2 [1 _ sftc^/i— -2 



(48) 



In the low-T regime, when either T -C min|e-|-| or T <C 7, the decoherence 
rate becomes linear in temperature: 
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-i(low) 



<r>-J\ 



4ir 



In 



1 - A(u> = 0) cosh 



OJ 

2T 



T . 2 f 2g 

— arctan 5— — - 

7r V4e 2 /7 2 - g 2 + 1 



(49) 



When both the coupling strength v and energy e are smaller than the tun- 
neling rate 7, the regime change occurs in the vicinity of T <~ 7 point (the 
upper curve on the insert in Fig. 7). If 7 is smaller than cither of the other 
two parameters: 7 <C ^4 = max{£,v/2}, then the high- and low-temperature 
asymptotes parametrically mismatch. In this case there is a crossover between 
the two regimes which is exponentially fast and occurs in the logarithmically 
narrow interval of temperatures: 



r 2 (x " over) (T)~r 2 (high) e - A / T t 



In A/ 7 



<T<A. 



(50) 



The temperature dependence of the decoherence rate is shown in Fig. 7 for 
the case of different couplings and fluctuator energies. We should stress that 




Fig. 7. Dependence of the decoherence rate on temperature, Eq.(47), for strong and 
weak coupling. The main picture shows a fluctuator with v e {uj) centered at e = 3, 
the insert shows e = 0. T, F2, e are measured in units of 7/2. 



although a linear in T behaviour similar to that in the Eq.(49) would also 
follow from the spin-boson models with the ohmic spectral function, as briefly 
described in Section 2, only a full quantum mechanical treatment can result 
in a non-trivial T dependence depicted in Fig. 7. 

Let us also stress that, although that at high temperatures decoherence 
is always stronger for a strongly coupled fluctuator, this does not have to 
be the case as temperature is lowered, since in the crossover region the g 3> 
1 fluctuators undergo stronger suppression than their j « 1 counterparts. 
Such a non-monotonic dependence of the decoherence rate J2 on the coupling 
strength v outside the classical region of high T deserves special attention. 
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When we start with a small value of v, v < 7, the peak of A(u>) grows as v 
increases which leads to higher decoherence rate r 2 . This proceeds on until the 
v = 7 point is reached when the peak stops growing but splits into two instead. 
One of those peaks then starts moving towards the origin which increases the 
overlap between the temperature function and A(u>), so F 2 keeps on growing. 
The maximum possible overlap is achieved when the incoming Lorentzian is 
centered at uj = which happens at the v — y/ 4e 2 + 7 2 point. If the coupling 
strength is increased beyond this point, the peak of A(u>) becomes de-tuned 
from the peak of the temperature function and the decoherence rate goes 
down. Since A(lo) plays the role of a spectral function, the positions of its 
maxima s± give us the energies of charge states of the system. Pushing them 
away from the Fermi level leads to freezing them out so that the impurity is no 
longer a fluctuator and thus does not contribute to the decoherence rate. Note 
that such a strong coupling is still not a good news for the qubit operation as 
the qubit and the impurity form together a four-level system which, although 
remains coherent, does not operate as intended. 

The above described behaviour of the decoherence rate as a function of 
coupling strength is depicted in Fig. 8. Note that the maximum of F 2 has a cusp 
which is not smeared by temperature (but would be smeared by including the 
a x part into the Hamiltonian) : Only at a rater high temperature (T <~ 100) 



Fig. 8. Non-monotonic dependence of the decoherence rate on the coupling strength 
at different temperatures for e — 3; here go = \fe 2 + 1. 

the decoherence rate practically saturates at its classical limit J^ hlg1 ^ =7/2. 
6 Relaxation in the FBC model 

Apart from decoherence, coupling of a qubit to the environment also leads to 
relaxation processes, i.e. to a decay of diagonal elements of the density matrix 
(23), as outlined in Section 2. In the model under consideration, Eq. (18), the 
presence of the cr^-part inevitably induces relaxation. Such a relaxation would 
also destroy coherence but the relaxation rate i~i is typically smaller than the 



T=100 



0.5 
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decoherence rate r 2 and therefore the time available for quantum operations 
is going to be determined by the latter. Experiments on the relaxation rate are 
of special interest as they can provide valuable information about the degrees 
of freedom in the environment to which qubit is coupled. 

It is convenient to diagonalize the qubit part of the Hamiltonian (18), 
transforming simultaneously the interaction part a z into a sum of f z and f x 
terms - the Pauli matrices in the new basis: 

H 



AE 

As it is well known [17], in Born-Markov rotating- wave approximation the 
decoherence rate A is given by: 

r 1= sin 2 6>S(AE), S(f2) = ([v(t)-(V),V(t')-(V)}) n , (52) 

where S(f2) is the spectral density of noise. (The appropriate derivation is 
outlined in Section 2). The symbol (. . .)q denotes the Fourier transform of 
the bath-averaged expression with respect to t — t 1 at the frequency i?. 

In the experiment of Astafiev et al [8] the authors determined the spectral 
density of noise in a wide window of frequencies, by changing AE and measur- 
ing the relaxation rate. The experimental data for S{fl) is presented in Fig. 9. 
Although the data points are somewhat scattered, the authors nevertheless 



AE 

— — t z + (cos 0t z + sin f x ) V + Hb ; where 

(5E c y + Ej, cos0 = ^|, sin0 = ^. (51) 




ta/2jr(GHz) 

Fig. 9. The experimental data by Astafiev et al [8]: the noise spectral density S 
deduced from measurements of A; open and closed symbols correspond to different 
samples while circles and triangles - to different measurement regimes; T~50mK. 



claim that they observed a linear behaviour of spectral density of noise with 
peaks at 7 and 23GHz. 

Below we will demonstrate that qualitatively the same behaviour is ex- 
pected from S(f2) in the FBC model at low temperatures. A contribution to 
the spectral density of noise from one impurity is given by 
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S(fl) = {#(t)d(t) - ($d) , d)(t')d{t') - M} ) = 

+ 00 

ttv 2 fl f 

= -^-coth— / &ujv e (u})v e (u} + fl) [n F (cu) - n F (uj + fl)} , (53) 

—00 

where the density of states v £ {uj) is given by (21), and np{u) is the Fermi 
distribution function. 

The frequencies accessed in the experiment [8] are in the window T <C fl < 
100T. Getting back to the rough estimates for the probability of having an 
effective impurity with energy smaller than temperature (the end of Section 
4), we see that even a scenario with \ei \ ^> 100T for all effective impurities is 
quite likely If this is the case then the condition T <C fl <C |e| holds and it is 
clear that the spectral density of noise S(fl) exhibits linear behaviour at such 
frequencies. As the thermal factor is practically equal to 1 within the range 
— fl < lo < and zero otherwise, the smearing of this step is of the order of 
temperature. Far to the right (assuming e to be positive) we have two peaks 
in DoS at e — fl and e. The main contribution to the integral (53) is coming 
from the tails of v e at the interval — fl < 10 < 0, rather than from the thermal 
function tail in the vicinity of DoS peaks which is exponentially suppressed. 
The latter contribution would only win if 7 itself was exponentially small, but 
the corresponding impurity is obviously non-effective. If 7 is not too small, 
we find 

2 

S(f2) = — ^(0)f2 T </2<e. (54) 

The summation of contributions (54) coming from all the effective impurities 
produces a linear in fl function, which can be a possible explanation for a 
linear (but very noisy) trend in S(£2) observed in [8]. 

Humps in Fig. 9 can be coming from less effective but more numerous 
impurities which makes the situation T <C \e\ < 100T possible. At the point 
fl <~ e the left one of the DoS double-Lorentzian overlaps with the thermal 
function thus producing a peak in S(fl). As fl further increases, the overlap 
between the left Lorentzian and the thermal function stays essentially the 
same but the region where it occurs drifts away to the left from the second 
Lorentzian which is stationary positioned at w = e. Qualitatively the picture 
of a linear behaviour of S(fl) followed by a peak at fl ~ s is the same for any 
relation between 7 and T, but the form of the peak is different. 

With a purely illustrative purpose we show in Fig. 10 a picture correspond- 
ing to the following scenario: all the tunneling rates are the same, 7, = 7, 
obeying 7 > T (which effectively allows one to put T to zero), and several 
impurities are uniformly distributed in some interval of energies. All these 
impurities contribute to the linear behaviour on the left of this interval [the 
grey area zoomed in the insert a)], while humps might be due to fluctuators 
with smaller energies weakly coupled to the qubit could not break the general 




Fig. 10. The noise spectral density S(fi) in arbitrary units at T = for the case of 
ten impurities with energies uniformly distributed between e = 60 and e = 240. The 
inserts zoom the bottom left corner (the grey area on the main figure for i < 40) 
when there is either a) no impurity in this interval or b) two impurities (si = 10, 
e — 20) with coupling constants 50 times weaker than those for 60 < e < 240. 



linear trend [the grey area zoomed in the insert b)]. The behaviour of S(f2) 
in the grey area qualitatively resembles the experimental data. 



7 Conclusions 

In these lectures, we have described some essential features of loss of coherence 
by a qubit coupled to the environment. We have first presented well known 
semiclassical arguments that relate both decoherence and relaxation to the 
environmental noise. Then we have shown that models with pure decoherence 
(when there is no relaxation in qubit states as that part of coupling to the en- 
vironment that leads to flipping the states is excluded) can be exactly solvable. 
As an example, we have treated in detail the model of fluctuating background 
charges [4-7] which is believed to describe one of the most important channels 
for decoherence for the charge Josephson junction qubit. Following our earlier 
treatment [7], we have shown that the decoherence rate saturates at 'high' 
temperatures [4, 5] while becoming linear in T at low temperatures and show- 
ing in all regimes a non-monotonic behaviour as a function of the coupling of 
the qubit to the fluctuating background charges. We have also considered, al- 
beit only perturbatively, the qubit relaxation by the background charges and 
demonstrated that the quasi-linear behaviour of the spectral density of noise 
deduced from the measurements of the relaxation rate can be qualitatively 
explained within this model in the low temperature regime. 

We thank B. L. Altshulcr, Y. M. Galpcrin, R. Fazio and A. Shnirman for 
useful comments. This work was supported by the EPSRC grant GR/R95432. 
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